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Introduction 

The term Yang-Baxter Equations (YBE) has been coined by Faddeev in 
the late 1970s to denote a principle of integrability, i.e. exact solvability, in 
a wide variety of fields in physics and mathematics. Since then it has be- 
come a common name for several classes of local equivalence transformations 
in statistical mechanics, quantum field theory, differential equations, knot 
theory, quantum groups, and other disciplines. We shall cover the various 
versions and their relationships, paying attention also to the early historical 
development. 

Electric networks 

The first such transformation came up as early as 1899 when the Brooklyn 
engineer Kennelly published a short paper, entitled the equivalence of tri- 
angles and three-pointed stars in conducting networks. This work gave the 





Figure 1: Star-triangle equation for impedances. 

definite answer to such questions as whether it is better to have the three 
coils in a dynamo — or three resistors in a network — arranged as a star or as 
a triangle, see figure d Using Kirchhoff's laws, the two situations in figure [T] 
can be shown to be equivalent provided 



ZiZi — Z2Z2 — ZsZs 

= Z1Z2 + 2^2-2^3 + Z-^Zi 

= Z{Z2Zs/(Z, + Z2 + Z3) 



[1] 

[2] 
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Here one has to take either POj or j2] as second hne of the equation, depending 
which direction the transformation is to go. The star-triangle transformation 
thus defined is also known under other names within the electric network 
theory literature as wye-delta (Y — A), upsilon-delta (T — A), or tau-pi 
(T — n) transformation. 

Spin models 

When Onsager wrote his monumental paper on the Ising model published 
in 1944, he made a brief remark on an obvious star-triangle transformation 
relating the model on the honeycomb lattice with the one on the triangular 
lattice. His details on this were first presented in Wannier's review article of 
1945. However, the star-triangle transformation played a much more crucial 
role in Onsager's reasoning, as it is also intimately connected with his elliptic 
function uniformizing parametrization. 

Furthermore, it implies the commutation of transfer matrices and spin- 
chain hamiltonians. Only in his Battelle lecture of 1970 did Onsager explain 
how he used this remarkable observation in his derivation of the formula 
for the spontaneous magnetization which he had announced as a conference 
remark in 1948 and of which the first complete derivation had been published 
by Yang in 1952 using a completely different method. 

Many other applications and generalizations have since appeared. Most 
generally, we can consider a system whose state variables — also called spins — 
take values from some suitable discrete or continuous sets. The interactions 
between spins a and b are given in terms of weight factors Wab and Wat, 
which are complex numbers in general, see figure El One quantity of special 
interest is the partition function — sum of the product of all weight factors 
over all allowed spin values. The integrability of the model is expressed by 
the existence of spectral variables — rapidities p, q,r, . . . — that live on oriented 
lines, two of which cross between a and b as indicated by the dotted lines in 
figure 121 Arrows from a to 6 are added to keep track of the ordering of a and 
b in case the weights are chiral (not symmetric). 

In Onsager's special Ising model case the spins take values a,b,c, . . . = 
±1 and the weight factors are the usual real positive Boltzmann weights 
depending on the product ab = ±1, uniformizing variable p — q and elliptic 
modulus k. In the integrable chiral Potts model the weights depend on a — b 
mod A^, with a,b = 1, . . . ,N, whereas the rapidities p and q are living in 
general on a higher-genus curve. 

When the weights are asymmetric in the spins, there are two sets of star- 
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Figure 2: Spin model weights Wab{p,g) and Wab{p,(l)- 

triangle equations which can be expressed both pictorially — see figure 0] — and 
algebraically: 

Wcd{p, q)Wdb{q, r)Wda{p, r) 

d 

= Rip, q, r) Wbaip, q)WUq, r)W,b{p, r) [3] 

R{p, q, r) Wah{p, qWac{q, r)Wbc{p, r) 

= ^^-(P' ^Wbd{q, r)Wad{p, r) [4] 

d 

Note that eqs. [H] and 0] differ from each other by the transposition of both 
spin variables in all six weight factors. In general there may also appear 
scalar factors R{p,q,r) and R{p,q,r), which can often be eliminated by a 
suitable renormalization of the weights. If a, b, and c take values in the same 
set, we can sum over a = b = c showing that R = R in that case. 

The Kennelly star-triangle equation P, [2j can be recovered as a special 
limit of a spin model where the states are continuous variables. 

Knot theory and braid group 

A seemingly totally different situation occurs in the theory of knots, links, 
tangles, and braids. In 1926 Reidemeister showed that only three types of 
moves suffice to show the equivalence between two different configurations, 
see figurelH Moves of type I — removing simple loops — do not apply to braids. 
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Figure 3: Star-triangle equation. 



Moves of type II, for which one strand crosses twice over another strand, can 
be reformulated for braids, namely that an over-crossing is the inverse of 
an under-crossing. The Reidemeister move of type III is a precursor of the 
more general Yang-Baxter moves and can be represented also by the defining 
relations of Artin's braid group. Let Ri,i+i be the operator representing the 
situation in which the strand in position i crosses over the one in position 
i + 1. Then a braid can be represented by a product of Rjj+i's and their 
inverses, provided 

Ri,i+lR«+l,j+2Ri,i+l — Ri+l,i+2R«,i+lRi+l,i+2 [5] 

and 

[R,,,+i,R,-,-+i]=0, if \i-j\>2 [6] 

and similar relations in which Ri,i+i and/or Ri+i_j+2 are replaced by their 
inverses. 

Factorizable S-matrices and Bethe Ansatz 

In the early 1960s Licb and Liniger solved the one-dimensional Bosc gas with 
delta-function interaction using the Bethe Ansatz. Yang and McGuire then 




Figure 4: Reidemeister moves of type I, II, and III. 

tried to generalize this result to systems with internal degrees of freedom 
and to fermions. This led to the discovery of the condition for factorizable 
S-matrices by McGuire in 1964, represented pictorially by figure El Here the 
worldlines of the particles are given. Upon collisions the particles can only 
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Figure 5: Vertex model Yang-Baxter equation. 

exchange their rapidities p, q, r, so that there is no dispersion. Also indicated 
are the internal degrees of freedom in Greek letters. In other words, the three- 
body S-matrix can be factorized in terms of two-body contributions and the 
order of the collisions does not affect the final outcome. McGuire also realized 
that this condition is all you need for the consistency of factoring the ra-body 
iJ-matrix in terms of two-body 5* matrices. The consistency condition is 
obviously related to the Reidemeister move of type III in figure |3] 
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Yang succeeded in solving the spin-i fermionic model using a nested Bethe 
Ansatz, utilizing a generalization of Artin's braid relations [S] and 

Ri,i+i{p - q)Ri+i,i+2{p - r)Rj,i+i(g - r) 

= Ri+i,i+2(g - r)Ri^i+i{p - r)Rj+i,i+2(p - q) [7] 

He submitted his findings in two short papers in 1967. The R operators 
in eq. [Zj — a notation introduced later by the Leningrad school — depend on 
differences of two momenta or two relativistic rapidities. Sutherland solved 
the general spin case using repeated nested Bethe Ansatze, while Lieb and 
Wu used Yang's work to solve the one- dimensional Hubbard model. 

Vertex models 

Since Lieb's solution of the ice model by a Bethe Ansatz there have been 
many developments on vertex models, in which the state variables live on 
line segments and a weight factors is assigned to a vertex where four line 
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Figure 6: Vertex model weight uj'^^{p,q), mixed model weight g) 
and IRF model weight w'^Kp, q). 

segments with the four states a, yU, A, jd on them meet, see figure El 

Baxter solved the eight- vertex model in 1971, using a method based on 
commuting transfer matrices, starting from a solution of what he then called 
the generalized star-triangle equation, but what is now commonly called the 
Yang-Baxter equation (YBE): 

5Z 5Z 5Z a" ^) '7" ^) 

a" 13" 7" 
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a" 13" 7" 

This equation is represented graphically in figure El From it one can also 
derive a sufficient condition for the commutation of transfer matrices and 
spin-chain Hamiltonians, generalizing work of McCoy and Wu who had earlier 
initiated the search by showing that the general six-vertex model transfer 
matrix commutes with a Heisenberg spin-chain Hamiltonian. To be more 
precise, Baxter found that if cu^^ = 5^5^ for some choice of p and q, some 
spin-chain Hamiltonians could be derived as logarithmic derivatives of the 
transfer matrix. 

Interaction-Round-a-Face model 

Baxter introduced another language, namely that of the "IRF-model" or 
"interaction-round-a-face" model, which he introduced in connection with 
his solution of the hard-hexagon model. This formulation is convenient when 
studying one-point functions using the corner-transfer-matrix method. Now 




Figure 7: IRF model Yang-Baxter equation. 

the integrability condition can be represented graphically as in figure [7| or 
algebraically as 

d 

= Y.'^d'iiP,<lX'aiQ,r)wfAP,r) [9] 

d' 

The spins live on faces enclosed by rapidity lines and the weights w'^Kp, q) 
are assigned as in figure El 
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Baxter discovered a new principle based on eqs. |Hland|ni which he called 
Z-invariance as it expresses an invariance of the partition function Z under 
moves of rapidity lines. This also implies that typical one-point functions 
are independent of the values of the rapidities, while two-point functions can 
only depend on the values of the rapidities of rapidity lines crossing between 
the two spins considered. Many recent results on correlation functions in 
integrable models depend on this observation of Baxter. 

IRF-vertex model 

In figureini we have also defined mixed IRF-vertex model weights W'^^\'^l{p, q). 
(We could put further state variables on the vertices, but then the natural 
thing to do is to introduce new effective weights summing over the states at 
each vertex.) With the choice made a more general Yang-Baxter equation 
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Figure 8: General Yang-Baxter equation, 
can be represented as in figure |S1 or by 

a" 13" 7" d 

x<4lSc.(g,r)<:^fe,:'(p,r) 

a" 13" 7" d' 

><Wt'^V:iQ,r)Wf;,\::^,ip,r) [10] 
Quantum Inverse Scattering Method 

The Leningrad school of Faddeev incorporated the methods of Baxter and 
Yang in their so-called Quantum Inverse Scattering Method (QISM), coining 
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the term Quantum Yang-Baxter Equations (QYBE) for the equations [S]. 
If special hmiting values of p and q can be found, say as h 0, such that 
5^5(^+0(/i), one can reduce |H] to the Classical Yang-Baxter Equations 



UJ 



A/3 



(GYBE) by expanding up to the first nontrivial order in expansion variable 
fi. These determine the integrability of certain models of classical mechanics 
by the inverse scattering method and the existence of Lax pairs. 

Checkerboard generalizations 

Star-triangle equations [2] and 0] imply that there are further generalizations 
of the Yang-Baxter equations, namely those for which the faces enclosed by 
the rapidity lines are alternatingly colored black and white in a checkerboard 
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Figure 9: Checkerboard versions of the weights. 

pattern. We can then introduce either vertex model weights u!'^^{p, q) and 
uj^'IIp, q), or IRF- vertex model weights W^'^^\'^l{p, q) and W^f^\'^l{p, q), or IRF 
model weights w'^Kp, q) and w'^{p, q), see figure IHl 

The black faces are those where the spins of the spin model with weights 
defined in figure |21 live; the white faces are to be considered empty in figures 
121 and El (or, equivalently, they can be assumed to host trivial spins that take 
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on only a single value). Clearly, the IRF- vertex model description contains 
all the other versions. 

Checkerboard vertex model 

First we consider the checkerboard vertex model with weights u^^{p, q) and 
uJ^^{p, q) as assigned in figure IHl The YBE |H] then generalizes to two sets of 
equations 

5^ 5^ 5^ f " (P> ^iy' 7 (P' 

a" /3" 7" 

{P,(l) 

a" f3" 7" 

x^:^,"{Q,r)u;];^;{p,r) [11] 

a" 13" 7" 
a" 13" 7" 

x^r7"(?,o<":(p,r) [12] 

where scalar factors R and R have been added as in jH] and jlj. These 
equations are represented graphically by figure UHl 

Checkerboard IRF model 

The checkerboard IRF version of the YBE |H] becomes 
Y,<'t'iP^<l)^tA<l,r)wt,::{p,r) 

d 

= R{p, q,r)J2 (P, q) <'i (Q, r) wf,, (p, r ) [13] 
d> 

Rip, ^iv (P, l) c' (p, r ) 

d 

= ^Va (P, l) <'i r) nfi\, {p, r) [14] 

d' 

again with scalar factors R and R added as in [H] and jlj. These equations 
can now be represented graphically as in figure ITTl Note that these equations 
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Figure 10: Checkerboard vertex model Yang-Baxter equation. 

reduce to eqs. [Hj and |lj if the spins on the white faces are allowed to take 
only one value, which means that they can be ignored. 

Checkerboard IRF-vertex model 

Finally, the most general case is represented by the checkerboard IRF-vertex 
model, with weights defined in figure IHl For this case the YBE are given by 




a" 13" 7" d 



a" 13" 7" d' 



[15] 



a" 13" 7" d 



x<4'il'(^'0<:a':'(p,r) 



11 



Figure 11: Checkerboard IRF model Yang-Baxter equation. 

a" (3" 7" d' 

xW^rfl^'l'(?,r)W^;qfi',(p,r) [16] 
with its graphical representation in figure El 
Formal equivalence of languages 
The square weight 

Combining four weights of a checkerboard model in a square, as is done with 
four spin model weights in figure ^1 we find a regular vertex model weight 
with rapidities that are now pairs of the original ones. In formula, this gives 

Wa^,{pl, qi)W^p{pi, q2)Wo,\{p2, gi)W^A/3(P2, q2) 

= ^af,iPi:P2;quq2) [17] 

From any solution of [Sj and ^ we can construct a solution of YBE |S] this 
way. This has been used by Bazhanov and Stroganov to relate the integrable 
chiral Potts model with a cyclic representation of the six vertex model. 
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Figure 12: Checkerboard Yang-Baxter equation. 
Map to checkerboard vertex model 

The checkerboard IRF-vertex model formulation contains all other versions 
mentioned above as special cases. However, collecting the state variables in 
triples, we can immediately translate it to a vertex model version, writing 

^liiP, Q) = WiXtiP, q), ZJ^ip, q) = WiXliP^ q) 

r A = (rf,A,c), /3=(6,/3,c) 
1^ d = (a, a, c?), fi = (a, /i, h) 

^&iiP^ 9) = ^lliP, = otherwise 

In eq. ^H] we have set all vertex model weights zero that are inconsistent 
with IRF-vertex configurations. Clearly, the translation of IRF models and 
spin models to vertex models can be done similarly. 

Map to spin model 

We can furthermore translate each vertex model with weights assigned as 
in figures IHl or El into a spin model with weights as in figure |21 by defining 

13 



[18] 
[19] 
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Figure 13: Square weight as vertex weight. 



suitable spins in the black faces, after checkerboard coloring. Each spin 
is then defined to be the ordered set of states on the line segments of the 
vertex model, a = [ai, a2, ■ ■ ■), ordering the line segments counterclockwise 
starting at, say, 12 o'clock. We can then identify uj^^{p,q) = Wa,b{p,q), 
uJ^^{p,q) — Wa^biP^Q.)- This is surely not very economical, as many of the 
weights will be equal, but it helps show that all different versions of the 
checkerboard YBE are formally equivalent. 

Hence, we shall only use the vertex-model language in the following. It 
is fairly straightforward to convert to the other formulations. 

An sl(m|n) example 

One fundamental example is a Q-state model for which the rapidities have 
2(5 + 1 components, p = ■ ■ ■ ,PQ),q= {q-g, ■ ■ ■,qQ), etc., and the states 
on the line segments are arranged in strings of continuing conserved color. 
The vertex weights, for a, P, X, /i — 1, . . . ,Q, are given by 



P+\ Q-13 
q+aP-n 



[20] 



with 



^^o%{Po: qo) = smh[r] + ep{po - go)] 
^o^piPo, Qo) = Gpa sinh(po - qo), P 
^o7p{po, qo) = e^°-i°>'^^P-''^ sinh 77, p a 
, '^Oa^^iPo, qo) = 0, otherwise 



[21] 
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where TV is an arbitrary overall normalization factor and rj is a constant. 
Furthermore, Ep = ±1 for p = 1, . . . ,Q, with m of them equal +1 and n 
of them equal —1. The Gp^s are constants satisfying Gpa = l/Gap, which 
freedom is allowed because the number of p-a crossings minus the number of 
cr-p crossings is fixed by the states on the boundary only, i.e. the choice of 
a, a', /3, /?', 7, 7' in YBE |H] and figure El 

The solution [20], |23 has many applications. The case m = 0, n = 2 
leads to the general six-vertex model; the m = 0, n = n case produces the 
fundamental intertwiner of affine quantum group \JqSl{n), whereas the case 
m = 2, n = l corresponds to the supersymmetric one- dimensional t-J model. 

Operator formulations 

The R-matrix 

For a problem with rapidity lines, carrying rapidities pi, . . . ,Pn, we can 
introduce a set of matrices Rij{pi,pj), for 1 ^ i < j ^ N, with elements 

R^Ap^,p,)t■:t = ^^A(p^,pj) n [22] 

In terms of these, the YBE jS] can be rewritten in matrix form as 

Rjk{Pj,Pk)RikiPi,Pk) Ri j (Pi , Pj ) 

= Rij{Pi,Pj)Rik{Pi,Pk)Rjk{Pj,Pk) [23] 

where l^i<j<k^N. 
The R-matrix 

If we transpose the P indices Pi and Pj in eq. |22], we can define a set of 
matrices Rj,j+i(p, g) with elements 

R^,+i{p,q)t:t = ^^^::M^) n [24] 

k^i,i+l 

Using these, the YBE jH] can be rewritten in matrix form as 

(g^ r ) Ri+i,i+2 (p, r ) Ri,i+i {p,q) 

= Ri+i,i+2(p,g)Ri,i+i(p,r)Ri+i,i+2(g, r) [25] 

and 

[RM+i(P.g).R.,.+i(^5)] = 0, if \t-j\>2 [26] 
In this formulation it is clear that many solutions can be found "Baxterizing" 
Temperley-Lieb and Iwahori-Hecke algebras. 
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Classical Yang-Baxter Equations 

If we expand 

RijiPi, pj) = 1 + hXijipi, p,) + 0{h^) [27] 

in we get in second order in h the Classical Yang-Baxter Equation 
(CYBE) as the vanishing of a sum of three commutators, i.e. 

\>(.ij{Pi,Pj),Xik{pt,Pk)] + \^ij{.pt,Pj),y^jk{,Pj,Pk)] 

+ \>^ik{PhPk),^3k{Pj,Pk)] = [28] 
introduced by Belavin and Drinfel'd, where Xjj is called the classical r-matrix. 
Reflection Yang-Baxter Equations 

Cherednik and Sklyanin found a condition determining the solvability of 
systems with boundaries, the Reflection Yang-Baxter Equations (RYBE), 
see figure Upon collisions with a left or right wall the rapidity variable 




Figure 14: Reflection Yang-Baxter Equation. 

changes from p to p and back. In most examples, in which the rapidities 
are difference variables such that R{p, q) = R{p — q), one also has p = fJ' — p, 
with n some constant. The corresponding left boundary weights are K^{p^p) 
satisfying 

Ki 1) Ri2 (p, q) Ki (p, p) Ri2 (g, p) 

= Ri2(p,g)Kib,p)Ri2(g,p)Ki(g,g) [29] 

with Ki(p, p) defined by a direct product as in |21] appending unit matrices 
for positions i ^ 2, and a similar equation must hold for the right boundary. 
Most work has been done for vertex models, while Pearce and coworkers 
wrote several papers on the IRF-model version. 
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Higher dimensional generalizations 

In 1980 Zamolodchikov introduced a three-dimensional generalization of the 
YBE, the so-called Tetrahedron Equations (TE), and he found a special 
solution. Baxter then succeeded in proving that this solution satisfies all 
Tetrahedron Equations. Baxter and Bazhanov showed in 1992 that this 
solution can be seen as a special case of the sI(cxd) chiral Potts model. Several 
authors found further generalizations more recently. 

Inversion relations 

When p) oc 5^5^^, i.e. the weight decouples when the two rapidities are 

equal, one can derive the local inverse relation depicted in figure E[ which is 




Figure 15: Local inversion relation. 

a generalization of the Reidemeister move of type II in figure jU It is easily 
shown that C{q,p) = C{p,q). 

This local relation implies also a global inversion relation which can be 
found in many ways. The following heuristic way is the easiest: Consider the 
situation in figure Uni with closed p-rapidity lines and M closed g-rapidity 
lines. For M and large, we may expect the partition function of figure 
ITHI to factor asymptotically in top- and bottom-half contributions. If each 
line segment carries a state variable that can assume Q values, then the total 
partition function factors by repeated application of the relation in figure El 
into the contribution of M + circles. Therefore, 

Taking the thermodynamic limit, 

z{p,q)^ hm ZM,^(p,g)i/*^^ [31] 
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Figure 16: Heuristic derivation of inversion relation. 



one finds 



z{p,q)z{q,p) = C{q,p) 



[32] 



In many models eq. |o2], supplemented with some suitable symmetry and 
analyticity conditions, can be used to calculate the free energy per site. 

See also 

Bethe Ansatz, Quantum groups, Ising model, Potts model. Vertex model, IRF model, 
Quantum inverse scattering method, Lax pair. Knot theory, Temperley— Lieb algebra, 
Factorisable S-matrix 



Further Reading 

Au-Yang H and Perk JHH (1989) Onsager's star-triangle equation: Master 
key to integrability. Advanced Studies in Pure Mathematics 19: 57-94. 

Baxter RJ (1982) Exactly Solved Models in Statistical Mechanics. Academic 
Press, London. 

Behrend RE, Pearce PA, O'Brien DL (1996) Interaction-Round-a-Face Mod- 
els with Fixed Boundary Conditions: The ABF Fusion Hierarchy. J. 
Stat. Phys. 84: 1-48. 



18 



Gaudin M (1983) La Fonction d'Onde de Bethe. Masson, Paris. 

Jimbo M (ed) (1987) Yang-Baxter Equation in Integrable Systems. World 
Scientific, Singapore. 

Kennelly AE (1899) The equivalence of triangles and three-pointed stars in 
conducting networks. Electrical World and Engineer 34: 413-414. 

Korepin VE, Bogohubov NM, and Izcrgin AG (1993) Quantum Inverse Scat- 
tering Method and Correlation Functions. Gambridge Univ. Press. 

Kulish PP and Sklyanin EK (1981) Quantum spectral transform method. 
Recent developments. In: Hietarinta J and Montonen G (cds) Integrable 
quantum field theories, Lecture Notes in Physics, vol. 151, pp 61-119. 
Springer, Berlin. 

Lieb EH and Wu FY (1972) Two-Dimensional Ferroelectric Models. In: 
Domb G and Green MS (eds) Phase Transitions and Critical Phenom- 
ena, vol. 1, pp 331-490. Academic Press, London. 

Onsager L (1971) The Ising model in two dimensions. In: Mills RE, As- 
cher E and Jaffee RI (cds) Critical phenomena m alloys, magnets and 
superconductors, pp xix-xxiv, 3-12. McGraw-Hill, New York. 

Perk JHH (1989) Star-triangle equations, quantum Lax pairs, and higher 
genus curves. Proc. Symposia in Pure Math. 49(1): 341-354. 

Perk JHH and Schultz GL (1981) New families of commuting transfer ma- 
trices in g-state vertex models. Phys. Lett. A84: 407-410. 

Perk JHH and Wu FY (1986) Graphical approach to the nonintersecting 
string model: star-triangle equation, inversion relation, and exact so- 
lution. Physzca A138: 100-124. 

Reidemeister K (1926) Knoten und Gruppen. Abh. Math. Sem. Hamburg. 
Univ. 5: 7-23. Elementare Begriindung der Knotentheorie. ibid. 24-32. 

Yang GN (1967) Some exact results for the many-body problem in one 
dimension with repulsive delta-function interaction. Phys. Rev. Lett. 
19: 1312-1314. 



19 



Yang CN (1968) S Matrix for the One-Dimensional A^-Body Problem with 
Repulsive or Attractive (5-Function Interaction. Phys. Rev. 167: 1920- 
1923. 

Remark: most older work can be found from Jimbo's book above, while 
most more recent work can be found at www.arXiv.org. 



Appendix: Relation of resistor networks with Gaussian model 
and Potts model 

We make next a few remarks on the relationship of resistor networks with 
related Gaussian models and Potts models in certain limits. 

We start with a graph G with vertices j and edges We asso- 

ciate with each vertex j an electric potential (pj and with each edge 
a resistance Rj,j'. According to Ohm's law the current from vertex j to a 
neighboring vertex j' (i.e. j' e Nbj) is given by 

^3,3' 

whereas the power dissipated in the resistors is 

{3,3') ^'^ 

with a sum over all edges. Kirchhoff 's second law is implicit in this formula- 
tion. However, minimizing P over the potential 0j for all internal vertices j 
(i.e. j e Into), we find 

and Kirchhoff 's first law emerges. 

We can now define the associated Gaussian model on graph G, using P 
for the interaction energy. The partition function is 

--(n;!/>'^.)-(-^i:^). i-i 
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which depends on the "potentials" of the exterior vertices. In the zero- 
temperature hmit (3 ^ oo the free energy tends to minP, defining the equi- 
hbrium state of the resistor network. 

We can even go further and derive a star-triangle relation for the Gaussian 
model. Doing a single Gaussian integration over 0o at the internal point 
of the star with vertices j = 0, 1, 2, 3, we find 

with j + 1 = 1 for j = 3 and Rj,j+i defined by 

3 



V^, for j = 1,2, 3. [38] 



The front factor in the RHS of [SZl has been factorized this way using the 
product of [211 over j = 1,2,3. It is related to the scalar factors R{p,q,r) 
and R{p,q,r) in [H] and The continuous variables (pj correspond to the 
discrete state variables a, 6, c, d and the integral a/z^/tt / d0o to the sum 
there. 

In the zero-temperature limit (3 —>■ oo the integral in the LHS of j2ZI is 
dominated by the maximum of the integrand and the RHS is also dominated 
by the exponential. The star-triangle relation P, [2] emerges, when iden- 
tifying Rjfl = Zj, {j = 1,2,3), and i?i,2 = ^3, ^2,3 = ^1, ^3,1 = ^2 in 



Another mapping has been given by Fortuin and Kasteleyn, who have 
shown that any planar resistor network is related to an — limit of a 
corresponding iV-state Potts model [See, e.g.. Section IV C of F.Y. Wu, Rev. 
Mod. Phys. 54, 235-268 (1982)]. Given Potts interaction energy 

we make, for each edge r = {j,j'), the special combinations 

_ e^-^- - 1 _ _ e^^- - 1 _ 1 



21 



The last equation defines the corresponding couphng constant Jr of the dual 
Potts model. The limit — is defined via the six-vertex model and in 
this limit we may choose {3 = Vn and find 

CC ip J y C I i^f^ ^ CC If J y C I ^ ^4 1 j 

with c some positive constant. 

Writing the star-triangle equation jS] for the iV-state Potts model in terms 
of X and x with 

Wab{p,q) = l + VNx{p-q)5a,b. Wab{p.q) = l + ^x{p-q)5a,b. [42] 

we find five independent equations. Eliminating R{p,q,r), we arrive at 

x{p — q)x{p — r)x{q — r) = x{p — q) + x{p—r) + x{q — r) + \fN, 

and x{u)x{u) = 1, [43] 

with the well-known solution x{u) = sin('u)/ sin(6' — m) = l/x{u), where 6 = 
arccos(v^/2). In the limit N this solution becomes x{u) = tan(M), 
x{u) = cot{u), while, in view of [H], Eq. [121 reduces to [T] and j2] with 

Zi = cta.n{p—q), Z2 = ccot{p—r), Z3 = ctan(g— r), 

Zi = ccot{p — q), Z2 = cta,n{p—r), Z^ = ccot{q — r). [44] 

This is a rapidity parametrization of the solution of and |2] . 
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